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C ' Abstract 

We study the probability densities of finite-time or local Lyapunov exponents 

(LLEs) in low-dimensional chaotic systems. While the multifractal formal- 
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ism describes how these densities behave in the asymptotic or long-time limit, 

(N 

O 

o. 

Q^ Depending on the nature of the dynamical state, the distribution of local 

Lyapunov exponents has a characteristic shape. For intermittent dynamics, 

>v 

and at crises, dynamical correlations lead to distributions with stretched ex- 

O . 

ponential tails, while for fully-developed chaos the probability density has a 

cusp. Exact results are presented for the logistic map, x — ► 4rc(l — x). At 

intermittency the density is markedly asymmetric, while for 'typical' chaos, 

it is known that the central limit theorem obtains and a Gaussian density 

results. Local analysis provides information on the variation of predictability 

on dynamical attractors. These densities, which are used to characterize the 

nonuniform spatial organization on chaotic attractors are robust to noise and 

can therefore be measured from experimental data. 



there are significant finite-size corrections which are coordinate dependent. 



I. INTRODUCTION 

The statistics of distributions of Lyapunov exponents (sometimes called stretch- 
exponents) have been studied in a number of physical situations ranging from turbulent 
flows || to Hamiltonian dynamics (in many-particle systems || and conservative mappings 
0), and are related to generalized dimensions and entropies @||. 

Lyapunov exponents provide a quantitative characterization of dynamics: for a dynam- 
ical system in D-dimensions, the D Lyapunov exponents of an orbit measure the rate at 
which volume elements in the phase space expand or contract along the orbit, in the dif- 
ferent directions. A positive Lyapunov exponent (LE) signifies exponential divergence of 
trajectories in the given direction, and is associated with chaotic dynamics, while negative 
LEs are associated with stable motion, when nearby trajectories converge. Although the 
exponents are global or asymptotic quantities, it is often instructive to examine the dis- 
tribution of values that the LE may take locally, namely over finite-time segments along a 
given trajectory [§-[jJ. If the underlying attractor is nonuniform, on a chaotic trajectory the 
local LE can be negative within a finite time interval. Similarly, on a nonchaotic trajectory 
the local LE can take positive values over finite time intervals ||. These considerations are 
of additional relevance for Lyapunov exponents computed from experimental data ||, since 
these are, naturally, only finite-time exponents. 

In this paper, we study the statistics of finite-time or local LEs in low-dimensional 
dissipative dynamical systems. Local Lyapunov exponents, which are defined in Eq. (^) 
below, depend on initial conditions, unlike the asymptotic or global Lyapunov exponent. 
The distribution of values that the local exponents take depends, in a characteristic manner, 
on the nature of the dynamical state; our present focus is on characterizing the different 
distributions that obtain for different dynamical attractors. 

We find that the characteristic densities of Lyapunov exponents fall in distinct classes 
depending on the nature of the attractor. Although the LLE distributions are stationary, 
they depend on the time interval over which the finite-time LE is computed. For very 



short times, the distributions keep changing shape and are difficult to classify, while in the 
asymptotic limit, as the time interval — > oo, all distributions must eventually collapse to a 
5-function centered on the global LE. The manner in which this happens is usually described 
by the multifractal formalism ||, but here we address the important corrections to scaling 
that can obtain for finite-times. 

Some aspects of such distributions have been studied previously. For "typical" chaotic 
dynamics, when correlations die out exponentially rapidly, the central-limit theorem holds 
for a number of averaged quantities, including local Lyapunov exponents |||7|]. Thus, the 
density is a Gaussian function, whose width depends on the length of the time interval over 
which the exponents are computed. For intermittent systems on the other hand, it is well 
known that correlations die out very slowly; this can lead to a power-law scaling for several 



quantities such as the Lyapunov exponent or the diffusion constant [10 1. Benzi et al. also 



studied intermittency in more detail and observed a deviation from the normal distribution 



in quantities such as the fluctuations of the response function JTT]|. 

We extend the analysis of finite-time Lyapunov exponents to the particular case of fully 
developed chaos and intermittency, where dynamical correlations persist over long times. 
This leads to significant departures from the simple central limit behaviour, and the resulting 
distributions are quite distinct from the normal density, typically having exponential or 
stretched exponential tails. 

An example which can be solved exactly is the commonly studied logistic map at the 
Ulam point, namely x —> 4x(l — x), for which we obtain an analytic form of the probability 
density for all times. The same distribution occurs for all parameter values where there is a 
boundary crisis, and thus appears to be quite general. 

We also treat the case of intermittent dynamics in some detail. In all instances of 
intermittency, the dynamics switches between two or more distinct types of behaviour. The 
distribution of LLEs that arise in such a situation can be shown to have components arising 
separately from these individual behaviours. 



Our main results are presented in Section II, where we discuss the different types of 
distributions for finite-time LEs. Our studies have been mainly on simple dynamical systems 
such as the logistic mapping, but the results we obtain appear to be more generally valid: 
these distributions can be seen a variety of systems (both mappings and flows). This is 
followed by a summary and conclusions in Section III. 

II. CHARACTERISTIC DENSITIES OF FINITE TIME EXPONENTS 



For generality, consider a .D-dimensional discrete nonlinear system |T2 



X n+1 = F w (X n ) (1) 

where X G R and {a} is a set of parameters. There are D Lyapunov exponents which are 
defined by considering a set of orthonormal D-dimensional vectors e m ,m = 1, . . . , D, and 
examining their evolution under the effect of the tangent mapping which is determined by 
the Jacobian of F, namely JF(X). Defining 

ef =(JF w (X J _ 1 ),ef_ 1 ), (2) 

the Lyapunov exponents are 

1 N 
A m = lim — yinlle"!, m = l,2,...D. (3) 

* t A 7" A J 11/11/ ? 1 \ J 
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The vectors e m are re-ortho normalized along the trajectory, and the subscript j here refers 
to the time. Stretch exponents are the logarithms of the ratios by which the vectors expand 
(or contract) along the D directions, 

Ar(j)=ln||e™||, (4) 

and this helps to define the m-th finite time Lyapunov exponent in a time interval of length 

iV as 

1 N 
ASf = T7EH|e?ll (5) 

JV 3=1 



In the remainder of this article only m = 1, the largest Lyapunov exponent, is considered, 
so the superscript index will be omitted henceforth in order to simplify notation. Local 
Lyapunov exponents are calculated along a trajectory that is divided into segments of length 
N. By Xn^) is meant the iV-step Lyapunov exponent calculated from the ith segment, and 
clearly 

i kN 

*N(k) = - £ Mj). (6) 

iV j=N{k-l)+l 

The asymptotic Lyapunov exponent A = Aqo does not depend (with probability 1) on initial 
conditions but A^ does. The probability density of local Lyapunov exponents, which is a 
stationary quantity, is defined as 

P a (X, N)d\ = Probability that Xn takes a value between A and A + dX. (7) 

If the stretch exponents of a system are considered to be random variables since the 
dynamics is chaotic, then the finite-time LEs should obey the central limit theorem, and 
the distribution function can be written in the general form 

P ( A - Af )~ [2,Afy]'/» B[p - Af * (A) - (8) 

where $(A) is a convex function with minimum at A = A. Expanding $ to second order 
gives the Gaussian density 

P Q (A,iV)~exp-(A-A) 2 /2a, (9) 

with variance a 2 oc 1/iV. Indeed this argument has been used very effectively to analyze finite 
time LEs on so-called "typical" chaotic attractors [|5|]F|,|10| as shown in Fig. 1(a), where the 
Gaussian nature of the density is evident. However, for other dynamical attractors, P(X, N) 
can be quite different, as shown in Figs. l(b)-(d). Indeed, the departure of $(A) from a 
polynomial function with quadratic maximum has been used to characterize the state and 
study the persistence of correlations [^, |T3i[T4| . 

Regardless of the behaviour of P(X, N) for small N, the local LEs eventually converge 
to the global exponent, \im.M-*oo P a (X, N) — > 5 (A — A). However for sufficient large A^ (but 



still far from the limit iV — > oo) the characteristic distributions depend on the details of the 
dynamics, so that the approach to the limit is distinctive. In the following subsections we 
discuss the particular cases of fully-developed chaos and intermittency in detail. 

A. Fully Developed Chaos 

The case of fully developed chaos in a system such as the logistic mapping, 

x n+1 = ax n (l - x n ) (10) 

at a = 4 can be analyzed in detail since the invariant density is known exactly. The 
Lyapunov exponent for this system is A = In 2, and the invariant density, which can be 
obtained by solving the appropriate Frobenius-Perron equation is [|l5|Jl6| 



p{x) = 1/k\Jx{1 -x),x e [0,1]. (11) 

The one-step Lyapunov exponent, 

Ai(i) = In |a(l- 2^)| (12) 

itself obeys the mapping 

Ai(i + 1) =ln|[exp{2Ai(i)}-a 2 + 2a]/2| (13) 

P(X, 1) is merely the invariant density for this mapping, and by using Eq. (0), one obtains 

P(A,1) = 2 ^P(^-^ 4 ) -oo<A<ln4. (14) 

7T A /l-exp[2(A-ln4)] 



P(X,N) can be calculated recursively since Xn is known in terms of Ai through Eq. (j^j. 
Since Ai is known in terms of Xi [Eq. (0)], it is possible to reexpress Eq. (|6|) as 

exp(iVA 7 v) - G(x) = 0, (15) 

where G is a polynomial function of x of order 2^-1. From this it directly follows that 



P(A,iV) = iVexp(iVA)£-^L ( 1 6 ) 

roots 1^ V X /I 

the sum being over all real roots, x(\), of Eq. fll5|) at given A G [—00, In 4]. Since the 
polynomial G(x) is of odd order, there is always one real root. For sufficiently small A all 
roots are real; with increasing A they leave the real axis in pairs, each such point giving rise 
to a singularity since the derivative dG(x)/dx vanishes there. 

In principle, Eq. (|i~6|) provides an exact solution for the invariant density of finite-time 
Lyapunov exponents for the logistic mapping. Furthermore, the same technique can be 
applied to obtain invariant densities for other systems if p(x) is known, though the form 
given above, namely Eq. (|T6|), is not particularly transparent. 

For this mapping, though, based on the analysis for small N, we conjecture the following 
asymptotic expression for the probability density of the iV-step exponent, 

«/* „ rt N exp(-JV|A-A|) , 

PA, N) = — ^ L- ^— — - 00 < A < In 4. (17) 

v ' ; n [1 - exp(-2iV|A - A|)]V2 - - v ; 

The main feature to be noted here is that there is a cusp at A = A, which does not vanish 
even as A" — > 00 although the range of A where it is significant decreases sharply with 
N. Outside this range, the function behaves essentially like an exponential. Eq. ( |17D thus 
provides finite-size corrections to the expressions derived earlier by Grassberger, Badii and 
Politi (see Eq. 4.9 in Ref. §). 

Both the results, namely Eq. ([UJ) or Eq. ( |iT|) can be verified numerically and can be 
shown to hold to a high level of accuracy. Shown in Fig. 2(a) are the (numerical) experimental 
distributions for the case of A^ = 14, which very closely matches both the exact (implicit) 
distribution, Eq. fll6|) , as well as the asymptotic result, Eq. [L7| [Fig. 2(b)]. Evaluation of the 
former expression requires the determination of the roots of the corresponding polynomial 
(we use the Newton-Raphson procedure). Shown in Fig. 3 is the solution for the case of 
N = 3, when the polynomial G(x) is of order 7. The divergences in the distribution P(X, 3) 
occur as pairs of real roots merge. For the specific case of the logistic mapping, as A" increases 
the largest number of singularities accumulate near A = In 2. The asymptotic form, Eq. ([Tj 



also gets progressively more accurate with increasing N [see Fig. 2(b)]. 

This form of the density is not restricted to the logistic mapping at a = 4, but is also seen 
in a number of mappings which have fully developed chaos. Even in the logistic mapping, 



it occurs at all parameter values corresponding to widening crises [0 , when the attractor 
is a rescaled image of the attractor at a = 4. 

We have also examined the dependence of the variance of the distribution on N. Typi- 
cally, a 2 oc 1/iV 7 for these distributions since they narrow with increasing N, going, in the 
limit of N — > oo, to a 5-function. For the case of a Gaussian density, 7=1, while for 
the exponential density the variance decreases more rapidly, and 7= 2. Our results for the 
variance are shown in Fig. 4. 

B. Intermittency 

Intermittent dynamics is characterized by a long-range temporal persistence of correla- 



tions Rig] , evidenced for example, by the existence of power-law dependence of a number of 
quantities on the parameters. The question of the distribution of local Lyapunov exponent 
for intermittent chaos has been explored previously by Benzi et al. [|ll|], who showed that 



there are significant departures from the Gaussian distribution. 

For the case of intermittency, the characteristic density of local Lyapunov exponent 
appears to be a combination of a normal density and a stretched exponential tail, an example 
of which is shown in Fig. 1(c). For A < A*, the density is a Gaussian, while above A*, the 
dependence is 

P(X, N) « exp[-N s (X - A*)] A > A,. (18) 

Since the exponent S < 1, the exponential tail decays extremely slowly with N. At A = A*, 
there is a crossover between the Gaussian and the stretched exponential density, which 
results in a completely asymmetric density about the mean. Similar distributions arise for 
all intermittent dynamical states, including the case of nonchaotic dynamics [|H| . 

8 



One way of understanding the above (phenomenological) expression for the density is 
to note that in all intermittent dynamics, the motion switches between two types of states. 
For each of these different dynamical states, the local Lyapunov exponents have a Gaussian 
distribution centered at different values of A, and with different amplitudes, and the stretched 
exponential behaviour interpolates between them. 

The example for which data is presented here in Fig. 1(c) and in Fig. 5 is the Type-I 
intermittency near the tangent bifurcation in the logistic mapping. This dynamics can be 
naturally separated into laminar regions (that stay close to the incipient period-3 orbit) and 
chaotic bursts. Finite-time Lyapunov exponents can be separately computed for trajectories 
that stay entirely in the laminar phase and entirely in the chaotic phase: these give the 
normal densities shown by dotted lines in Fig. 5. Trajectory segments that visit both the 
components of the intermittent motion contribute to the stretched exponential tail; with 
increasing N, the purely chaotic component is more difficult to identify, since there are 
fewer segments that are of duration longer than N. The stretched exponential tail thus 
decreases with increasing N, and the distribution eventually collapses to a delta-function. 

This behavior is generic at all intermittencies. For the mapping 

which has been extensively studied within the thermodynamic formalism in the context of 
Type-I intermittency |2(J , we observe a similar decomposition of the overall density of LLEs 



to a superposition of two independent Gaussians with stretched exponential interpolation 
between the two. Furthermore, we have also examined a number of higher dimensional maps 
and flows, and find that at all intermittent dynamics, including the cases of forced systems 
or of nonchaotic dynamics |Ti9|j2T| , non-Gaussian stretched exponential tails are seen. An 
example shown in Fig. 6 is the density P(\, 2048) for the largest nonzero Lyapunov exponent 
in the Lorenz system (see the figure caption for a definition of the dynamical system) with 
the parameters chosen to correspond to intermittent dynamics. 

The case of crisis-induced intermittency @ (just beyond widening crises, for instance) is 

9 



similar, with the two independent densities being, respectively, the exponential cusp, namely 
Eq. (|17D for the component for the pre-crisis chaotic attractor, and a Gaussian density which 
corresponds to the widened chaotic attractor. The density shown in Fig. 1(d) for the period- 
five widening crisis in the logistic mapping can be analysed in a manner very similar to the 
case illustrated in Fig. 5. 

The variance for these distributions decrease somewhat faster that for the Gaussian, 
namely a 2 oc 1/A 7 with 7 = 25 > 1. With increasing N, the exponent however changes, 
eventually reaching the Gaussian limit, 7 « 1 (see Fig. 5). 

III. SUMMARY AND DISCUSSION 

A major motivation for the present work has been the realization that local Lyapunov 
exponent distributions have characteristic forms depending on the nature of the attractors, 
and that these provide an additional and important dynamical characterization of the chaotic 
state of a system. 

We have mainly focused on the cases of fully developed chaos, crises and intermittencies — 
namely those attractors for which the density shows a marked departure from the simple 
Gaussian form. This is indicative of significant finite-size corrections to the multifractal 
formalism ||22| . All these cases show exponential tails albeit for different reasons |23|]. For 
fully developed chaos in the logistic mapping we obtain the (in principle) exact expression 
for the probability density as well as an asymptotic approximate form. These densities are 
seen at all parameter values corresponding to interior crises, and thus are quite common. 

The case of intermittency was analysed in detail, and the density seen there was shown 
to arise from individual densities corresponding to the different components of the motion 
(laminar and chaotic, say) with the stretched exponential tail corresponding to interpolation 
between these two. 

We have verified the generality of our results for a number of other dynamical systems, 
for example higher dimensional mappings such as the Henon system or flows such as the 
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Duffing oscillator, the forced damped pendulum, the Lorenz equations etc. These densities 
are quickly attained, and are maintained even as N increases, and for fairly high levels of 



additive noise [24| . This is of particular relevance when analyzing experimental data. 

The logistic map was used here for illustration since exact results can be obtained for at 
least one parameter value, but the probability density for finite-time Lyapunov exponents 
in any system can be obtained via Eq. (|T6| ) so long as the invariant measure is known. In 
this regard, it is interesting to note that recently, Pingel et al. |26j have described a general 
inversion technique whereby a class of 1-d maps having a prescribed invariant density can 
be constructed. 

The characteristic forms for the density that we have described are found in a variety 
of systems, including those where the dynamics is not chaotic. When a system is forced 
quasiperiodically, chaotic attractors can be transformed to strange nonchaotic attractors 
||2 1|| . These attractors are fractal, but the largest Lyapunov exponent is zero or negative. 
Because of the spatial fractality, though, for short times the local Lyapunov exponents can 
be positive, and the distributions again fall into the classes that we have seen here for 



chaotic systems |27| . The phenomenon of high-stretch tails also appears to be very general. 
In recent work Calvo and Labastie |28j have examined the distribution of local Lyapunov 
exponentsin a conservative system, namely in a 19-particle cluster simulation study. They 
observe that when the dynamics is intermittent |29|| , the Lyapunov exponent density has a 
stretched exponential tail. 

Local analysis can be more revealing of the nature of the dynamics than global quantities, 
and this is an issue when predictability is of concern. For instance, in applications that aim 
to predict future behaviour based on time series data (for example in atmospheric sciences 
or in economics), atypical or extreme behaviour which contributes to the stretched tails is a 
serious bottleneck. Since the largest Lyapunov exponent can be extracted from time-series 
data though standard techniques, examination of the finite time distributions can give more 
insight into the dynamics than the extraction of a single exponent. 
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We conclude with a few general comments. Except for fully developed chaos when the 
distributions can be derived for all times, the present study does not examine the case of 
very short times when the distributions are atypical and it is not clear if they are stationary. 
The characteristic behaviour becomes apparent for times that are not too short, and persists 
thereafter. Thus, local Lyapunov exponent densities provide quantitative distinctions among 
different chaotic attractors. Exponential tails are characteristic of fully developed chaos and 
intermittency: whether these are related to analogous distributions that arise in turbulent 
flows ]30| or stretched-exponential tails in relaxation phenomena is an interesting open 



question. 
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Figure Captions 

Fig. 1 Characteristic probability densities that arise in chaotic dynamical systems, corre- 
sponding to (a) "typical" chaos, (b) fully-developed chaos, (c) intermittency, and (d) 
crisis-induced intermittency. Numerical results are for P(X, N) for the logistic map 
with (a) a = 3.7, (b) a = 4, and (c) a = 1 + v^-lO" 6 , and (d) a = 3.7 447 104 . The 
local Lyapunov exponents are calculated from N step segments with 10 6 iterations. 
Different values of N are chosen for the three cases, in a) the Gaussian nature becomes 
apparent at N around 10 or so, the cusp in b) is already evident for N around 8 and 
survives for all N, while the asymmetry in the distribution c) and d) appears to persist 
for all N. Note that for d), N should be a multiple of the periodicity of the window. 

Fig. 2 (a) Comparison, at a = 4, of the numerical results (dots) for the density P(\, 14), with 
the analytic expression, Eq. flilf) , (solid line), (b) Numerical results (dots) compared 
with the approximate density, Eq. fll7|) (solid line), for N = 100. 

Fig. 3 The real roots (solid line) and P(X, 3) (dotted line) for the logistic map with a = 4. The 
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singularities in the density occur when a pair of real roots become complex. P(X, 3) 
has been rescaled for clarity. 

Fig. 4 (a) The variance as a function of iV for "typical" chaos (□), fully-developed chaos 
(o), and crisis-induced intermittency (V). The exponents characterizing the decay 
are 1.12,1.95, and 1.51 respectively. For intermittent chaos (filled circles) there is a 
crossover: the exponents in different ranges of N going from exponential limit, 1.85 to 
the Gaussian limit 1.09 at large N. 



Fig. 5 Near the intermittent transition, at [a — 1 + \/(8) — 10 -6 ], where there are long 
crossovers, the two components of the density P(X, 300) (dotted lines) are compared 
with the total density (solid line). For this state, X p = —0.00015 and A* = 0.03. Note 
that the amplitudes of the individual densities have been appropriately scaled to depict 
clearly the manner in which they contribute to the total density. 

Fig. 6 For the Lorenz equations : x = a(y — x), y = x(r — z) — y,z — xy — bz, where a, 
b and r are parameters, the density P(X, 2048) near the intermittency transition at 
a = 10, b = 8/3 and r = 166.8 801 548. Here A is the largest nonzero Lyapunov 
exponent, and the integration step size is 0.02 natural units. 
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